O 
O 
(N 

u 

a* 
< 



O 

CO 



c3 



Correlated and uncorrelated nanoscale heterogeneities in Ll solid solutions and their 

signatures from local and extended probes 
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The phase coexistence of chemically ordered Llo and chemically disordered structures within 
binary alloys is investigated, using the NiMn system as an example. Theoretical and numerical 
predictions of the signatures one might expect in data from local and extended probes are presented, 
in an attempt to explain the presence of antiferromagnetism in NiMn when no Llo signatures appear 
in diffraction data. Two scenarios are considered, the first in which the tetragonal Llo structure 
and fee chemically disordered structure are distributed evenly into uncorrelated domains of specified 
average diameter. The diffraction limit, below which the two structures can only be distinguished 
using a local probe, is quantified with respect to the domain diameter by applying straightforward 
diffraction analysis. In the second scenario, domains with chemical ordering oriented in different 
directions are required to maintain their atomic coherence with each other. A numerical treatment 
is used to illustrate the long-range strain that results from elastic energy considerations, and the 
effects on the structure factor (extended probe) and pair distribution function (local probe) are 
investigated. 
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I. INTRODUCTION 

The M X F\- X intermetallics, where a metal, M, is com- 
bined with F, one of the three ferromagnetic transition 
metals iron, nickel, or cobalt, for x close to 0.5 or 0.75, 
provide some of the best known examples of the im- 
portance of chemical ordering as determinants of struc- 
ture:property relationships. When crystallographically 
disordered by being prepared as, e.g., thin films deposited 
on low temperature substrates, they are paramagnetic 
and cubic. However, when they are allowed to develop 
their most stable structures by preparing or annealing 
them at higher temperatures, the different elements or- 
der as homogeneous layers occupying the [100] planes 
that alternate with each other according to the overall 
stoichiometry of the material. They are then commonly 
referred to as Llo alloys because the different sizes of 
their constituent elements as well as the chemical order- 
ing results in a tetragonal distortion perpendicular to the 
layers. The proximity and periodicity of the atoms now 
cause the spins on the atoms of the FM metal to order 
first within a layer and next between the neighboring lay- 
ers so that the Llo alloys are ferro or antiferromagnets, 
with the magnetization also along the tetragonal axis. 
Of obvious interest are the intermediate regime between 
the fully disordered and ordered materials, the path con- 
necting them, and the relationship between the partial 
structural spin ordering. 

In the case of Ni^Mni-^ alloy films near the equiatomic 
composition x — 0.5, the magnetic and structural prop- 
erties have generated interest both from the basic and ap- 
plied viewpoints. Such films exhibit the aforementioned 
LIq tetragonal phase which comprises alternating atomic 



planes of Ni and Mn along the c-axis with a contraction 
perpendicular to these planes. In its chemically ordered 
state, NiMn is a strong antiferromagnet, which makes it 
a good candidate for use as the exchange coupling mate- 
rial to pin the reference ferromagnetic layer in spin valve 
magnetic read sensors i£ Chemically ordered NiMn is ca- 
pable of providing strong pinning fields, a high blocking 
temperature, and good corrosion resistance^*^*^ 

It has been shown that the exchange bias in NiMn 
alloys decreases with decreasing antiferromagnet layer 
thickness and that a critical thickness is needed to de- 
velop the exchange coupling. It has been suggested^ that 
such thin films (< 5 — 10 nm) are not chemically ordered 
and hence are not antiferromagnetic, which explains the 
small exchange coupling in these films. In grazing in- 
cidence X-Ray diffraction measurements up to Q ~ 5 
A -1 on NiMn films annealed at 280° C over this range 
of thicknesses, no traces of sites were found. In 15 — 25 
nm thick films, the chemically ordered Llo phase grows 
at the expense of the fee phase, increasing the exchange 
bias to 300 - 400 Oe. 

An intriguing aspect of these M X F\- X intermetallics is 
that the onset of magnetization occurs prior to the obser- 
vation of chemical ordering/tetragonal distortion. This 
behavior suggests that the spins can align in magnetic 
domains that are at or below the size of structural do- 
mains that would give signatures in diffraction, i.e., the 
diffraction limit. Local structure measurements on PtCo, 
PtFe, and NiMn have in fact demonstrated a correlation 
between this initial magnetization and ordering on local 
length-scales i^iiSi It is notable that a very high degree 
of local ordering was observed in these studies that sug- 
gests the materials are actually composed of incoherent 



nanoscale domains that are close to the homogeneous Llo 
structure. 

XAFS studies provide complementary information not 
only to diffraction data, which gives the long-range aver- 
age order, but also to that obtained from diffuse scatter- 
ing and pair distribution analysis, providing an alterna- 
tive means for probing the development of the nucleation 
sites for the Llo phaseiii*i£ Indeed, XAFS has been em- 
ployed at the Mn K-edge as a direct and element specific 
probe of the average local atomic structure within ^5 — 6 
A of the Mn site, in which a short-range order parameter 
is introduced to quantify and confirm the presence of Ll 
structures^ At low annealing temperatures, the short- 
range order parameter differs from an associated long- 
range order parameter obtained from diffraction data, 
providing a sufficient condition for nanoscale phase sepa- 
ration and heterogeneity by demonstrating the presence 
of a nondiffracting but locally ordered component within 
the alloy. In these films a nonzero coercivity is measured 
even while the associated diffraction patterns exhibit no 
signs of an Llo structure. 

When correlations between local distortions produce 
ordered distributions, they are generally observed in 
diffraction experiments as superlattice or satellite peaks 
(in addition to the Bragg peaks) >i» The situation, how- 
ever, is more complicated when domains below the limit 
of long-range order occur in a random or aperiodic 
distributionii£*i£ii£ In this case, the effects of heterogene- 
ity, defect or other chemical organization, and competi- 
tion between different phases on the diffraction patterns 
are not well understood. In Garcia- Adeva et al.mi for ex- 
ample, it has been shown that a second coexisting struc- 
ture affecting one-fifth or more of the atoms in various 
two dimensional lattices can have a minimal to negligible 
effect on the calculated structure factors when compared 
with their periodic counterparts. The signatures of dis- 
order that are present in the diffuse part of the scatter- 
ing are usually smooth and subtle, and it is often easy 
to construct several distinct disordered lattices that pro- 
duce essentially equivalent diffraction patterns. 

The indications of chemical ordering within very small 
domains of sizes that will average crystallographically to 
a random, cubic structure is not surprising. Although 
even a random distribution of atoms may produce a 
plethora of semi-ordered domains containing up to sev- 
eral dozen atoms, given the stability of the layered struc- 
tures in this system it is likely that when new atoms 
are deposited they add to existing structures in an or- 
dered way that is terminated when the growing layered 
nanocrystallites collide. An alternative description might 
be to consider a random solid solution where each atom 
is allowed to exchange its position with a neighbor to 
form the more stable ordered structure only once, result- 
ing in a large number of small domains that individu- 
ally exhibit a high degree of order internally. Because 
of the significance of this phenomenon whereby various 
properties have size thresholds for their expression in this 
and other systems, it is worth exploring their structural 



properties so as to understand the experimental signa- 
tures of nanoscale ordering in crystalline materials and 
the ramifications of ordered domains on stress and defect 
formation. 26 

Here, using the example of the NiMn system, we inves- 
tigate the effects that domain size and interface behavior 
have on conventional diffraction when nanoscale phase 
separation and heterogeneity either occur below such a 
diffraction limit, or are faced with the restriction of hav- 
ing coherent atoms across their respective interfaces in 
which long-range elastic strain prohibits their preferred 
structures from forming within the domains. We study 
the relevant length-scales in which chemical order can ex- 
ist in the form of small nanoscale Llo domains immersed 
in a chemically disordered matrix with an average, crys- 
tallographic structure, and quantify the domain sizes by 
which such distortions will not be apparent in the diffrac- 
tion data and yet still resolvable with local probes. In ad- 
dition, at this length-scale, we discuss the role that elas- 
tic stresses at the domain interfaces play in determining 
the relaxation field of the structure within the domains 
themselves, the energy considerations of such interfacial 
interactions, and the associated signatures such a strain 
has on the long-range (e.g. diffraction) and local (e.g. 
XAFS, pair distribution) probes. 

The analytical and numerical analysis that follows will 
make use of the structure factor S(Q) and pair distribu- 
tion function (PDF) G(r) of various lattices to illustrate 
the effects of nanoscale heterogeneity on extended and 
local probe data, as the two functions are well-defined 
and computationally tractable. The formalism used is 
identical to Peterson et aliU with comparisons to other 
definitions and nomenclature found in Keen 28 . 

Nanoscale heterogeneity and phase separation are not 
exclusive to the magnetic alloy considered in here. This is 
also the case, for example, with substitutional impurities 
in solid solutions. 19 If the impurity atoms attract or repel 
each other, heterogeneous nanoscale domains or texture 
will form in the solid. Likewise, these effects are seen in 
charge separation and charge ordering in colossal mag- 
netoresistance compounds and high T c superconductors, 
leading to the formation of stripesi 20 i 21 i 22 i 23 ' 24 i 25 



II. INTRINSIC NANOSCALE ORDERING AND 

THE FORMATION OF CORRELATED AND 

UNCORRELATED DOMAINS 

A random arrangement of atoms will exhibit chemical 
ordering on all length-scales, with more ordering appear- 
ing more often at smaller, nanoscale lengths. In addi- 
tion, a growth process might favor the formation of such 
domains beyond their intrinsic random existence, if the 
energetics permits. The growth of layers during atom de- 
position, or neighbor exchange on a surface, are examples 
mentioned before. 

In the case of the Llo ordering just described, one 
might define a region, or domain, of such ordering if 
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FIG. 1: A random arrangement of red and blue atoms will 
always display some ordering, depicted as rectangles in this 
figure. In certain cases, such ordering can trigger the onset of 
a function, like antiferromagnetism. This is illustrated with 
the use of an order parameter 0(x) that attains some critical 
value within the cores of ordered domains. 



more than 3/4 of the atoms satisfy the definition of Llo 
ordering, with alternating rows of distinct atoms. Fig- 
ure ^ shows a two dimensional example of this, in which 
equal numbers of red and blue atoms are distributed ran- 
domly. Rectangular regions are demarcated if their en- 
closed atoms show such ordering, and the parallel lines 
within the domains show the LIq orientation. Using 
these rules, it is seen that roughly 35 — 40% of the atoms 
are found in such layered-type domains. With this, one 
can define an order parameter 0(x) that has a value of 
zero in locations where no LIq ordering is present, a value 
of one within domains exhibiting LIq ordering, and a 
value of one-half at the interfaces. The top of the Fig. 
^ shows a plot of cf> across the middle horizontal row of 
atoms, in which the ordering is high (solid line), and a 
horizontal row from a lower portion of atoms, in which 
the ordering is low (dashed line), along with values of the 
average ordering avc and a hypothetical critical value 
0crit- When the ordering is above the critical value, the 
structure takes on a function, in this case antiferromag- 
netism. Recent calculations of local electronic density of 
states in NiMn do in fact show that the perturbation on 
the spin ordering from the inhomogeneity at the inter- 
face is spatially confined . 2 - The plot shows that indeed 
certain regions attain this value, despite the random ar- 



rangement of atoms. The majority of ordering will exist 
at the nanoscale, and hereafter we set out to determine 
the signatures that such ordering exhibit in extended and 
local probes. 

Our treatment of the fcc-Llo phase coexistence in 
NiMn alloys will mostly involve cases in which the two 
structures are present in equal amounts, and we consider 
the case of total LIq ordering for completeness. We have 
found that the degree of atomic correlation at the inter- 
faces of domains plays a critical role in determining the 
actual crystallographic structure within the domains. In 
the case where there is no correlation whatsoever between 
the orientation and placement of one domain with respect 
to another, then it is an excellent approximation to as- 
sume that the internal structure of the domain is iden- 
tical to that of an associated infinite crystal (that is, it 
is "ideal"). Small distortions at the interfaces contribute 
insignificant signatures to S(Q) and G{r). 

For small domain sizes, the structure within the do- 
mains can be close to ideal if one allows for dislocations 
at the interface between the domain and the matrix it 
is embedded within. Dislocations serve to localize the 
large elastic energy present in any system containing two 
or more phases with different lattice parameters between 
them. Without dislocations, the internal structures of 
the domains and the matrix would be far from ideal, ex- 
hibiting a continuous strain throughout such that the 
structure remains coherent at the interfaces. 

The top image in Fig. [5] illustrates how a disloca- 
tion at the interface can localize the strain such that the 
internal structures of LIq domains and the fee matrix 
could take on their ideal structures. The [111] plane of 
a NiMn alloy with both LIq and fee structures can be 
described as a triangular arrangement of atoms, with the 
LIq (tetragonal) structure in this plane having expanded 
atomic separations along rows of similar atoms and con- 
tracted atomic separations along rows of different atoms. 
The fee structure has a higher symmetry chemically dis- 
ordered phase in which the average bond length between 
any two neighbor atoms is constant. The lattice mis- 
match between the structures is accounted for locally by 
the presence of a dislocation, shown in the inset. Al- 
though the dislocation introduces a large distortion lo- 
cally in the lattice, distant atoms are able to take on 
their preferred arrangement, even at the interface of the 
LIq and fee domains. 

The bottom image in Fig. |5]depicts an arrangement of 
hexagonal domains with complete LIq ordering, in which 
no dislocations are allowed. The resulting stress on the 
atoms induces a strain that propagates throughout the 
domains, as depicted in the inset. The length-scale over 
which the strain extends tends to alter the unique tetrag- 
onal structures that are otherwise present in the case of 
the rectangular domains that allow for dislocations. We 
will see that such an arrangement of domains, and the 
accompanying strain, can actually produce a structure 
with an increased tetragonal distortion, albeit highly dis- 
ordered. 
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FIG. 2: Top: A [111] plane showing portions of three chem- 
ically ordered Llo domains separated by a chemically disor- 
dered fee phase. The lattice mismatch is accounted for locally 
by the presence of a dislocation (inset), which might allow 
the structure to attain the separate fee and tetragonal phases 
elsewhere. Bottom: complete Llo ordering within hexagonal 
domains can force the crystal into a highly strained state if 
no dislocations are allowed, thus diminishing the presence of 
either fee or tetragonal structures. 



It is useful to consider whether the presence of disloca- 
tions, which localize the large elastic energy around a core 
of atoms while allowing a minimum elastic energy config- 
uration elsewhere, is energetically favorable over the case 
of an identical distribution of domains without disloca- 
tions. The latter scenario tends to uniformly distribute 
the elastic energy throughout a large region near the in- 
terfaces. Consider a [111] plane of NiMn alloy in which an 
arrangement of two dimensional LIq rectangular domains 



is embedded within a chemically disordered fee structure, 
as depicted in the upper half of Fig. [21 The total elastic 
energy is minimized by varying the positions of the atoms 
(the exact interaction potentials between the atoms are 
detailed in the fourth section). In one case, dislocations 
are allowed at the interfaces. In the other case, the dis- 
tribution of domains is identical, but no dislocations are 
present and each row of atoms spans the entire lattice, 
thus maintaining coherence. 
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FIG. 3: The average energy per atom as a function of rect- 
angular domain size, characterized by the average number of 
atoms N per side. The lattice mismatch of the two struc- 
tures determines the domain size in which dislocations at the 
interfaces are favorable. 

Figure [3] shows the minimum elastic energy of uncor- 
rected (dislocations allowed, incoherent) and correlated 
(no dislocations, coherent) domains for the two dimen- 
sional case. The rectangular domain size is characterized 
by the average number of atoms N per side. Clearly, 
it becomes more favorable to have uncorrelated domains 
coexisting within the two dimensional structure as the 
average domain size increases. For N < 100, however, 
dislocations at the interfaces introduce more elastic en- 
ergy per atom than highly correlated structures in which 
the strain is uniformly distributed throughout the inte- 
rior of the domains. In this case, the lattice mismatch 
between the Ll and fee structures is too small to be 
accounted for by the addition or subtraction of one dis- 
location row. 

This analysis of a two dimensional system motivates 
one to consider both scenarios. Clearly, correlated do- 
main distributions are favorable when the domain size 
is small, and we will investigate their effects on S(Q) 
and G(r) for various domain sizes. In this case ener- 
getics plays an important role, and the atom positions 
(and the ensuing long-range strain) must be carefully ac- 
counted for with a robust numerical model. In the next 
section, however, we consider the case in which spherical 
domains of ideal LIq and fee structures coexist in a three 



dimensional structure. By incorporating the reasonable 
approximation that the domains are entirely uncorrelated 
with each other, through the presence of dislocations or 
otherwise, we quantitatively obtain their diffraction and 
PDF signatures with little analytical effort, and without 
any need for large-scale numerical treatments. 



III. THE STRUCTURE FACTOR OF THREE 

DIMENSIONAL UNCORRELATED FCC AND Ll 

SPHERICAL DOMAINS 

The radially averaged structure factor, S(Q), of a 
NiMn alloy with uncorrelated fee and LIq spherical do- 
mains can be derived analytically with very little approx- 
imation needed. Consider a NiMn alloy in which the fee 
domains have an internal structure defined by a unit cell 
with lattice parameter a = 3.65 A, and the Llo tetrago- 
nal domains have an internal structure defined by a unit 
cell with a = 3.72 A and c = 3.52 AM The atomic densi- 
ties of the two internal structures are equal within 0.2%. 

The structure factor can be obtained directly from the 
pair distribution function, G(r), by a Sine transform. 
We have previously shown that uncorrelated spherical 
domains with a particular crystal structure and aver- 
age diameter d generate a PDF that is simply that of 
the associated infinite crystal multiplied by an envelope 
function that is characterized by dw& As a result, S(Q) 
for the uncorrelated spherical domains is obtained from 
that of an infinite crystal by a related convolution, that 
broadens the Bragg peaks, giving them a half-width of 
approximately 



AQ = 



3.48 



(1) 



Consider the effect such a convolution has on the ideal 
Bragg peaks in the interval < Q < 5 A -1 . The con- 
clusions made within this interval can easily be extended 
to larger intervals. The ideal fee structure in the NiMn 
alloy has three Bragg peak positions at Q\ = 2^/3tt/o,, 
Q2 = Air/a, and Q 3 = Any/2 /a. The ideal Ll struc- 
ture has five Bragg peaks at Qi, Q2 — AQ^, Q2 + AQ^, 
Q3 — AQ3 , and Q3 + AQ3 , where the splitting widths 
from the second and third fee peaks are given by 
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in S(Q). As the average domain size decreases the broad- 
ening of the peaks introduces an overlap between neigh- 
boring peaks that eventually removes each Bragg peak 
maxima from the LIq structure (as its peaks are smaller 
in magnitude than those of the fee structure) . The peaks 
then become "shoulders" of their neighboring fee peaks, 
and are still distinguishable because of changes in the 
local curvature of the structure factor between the two 
peak locations (i.e., an inflection point still resides be- 
tween them). Eventually, however, the average domain 
diameter d becomes small enough that such curvature 
disappears altogether, and the LIq peak is no longer re- 
solvable. 

Such a "diffraction limit" for the domain diameter can 
be easily determined by equating the convolution half- 
width given by Eq. ^with each of the splitting widths in 
Eq. [5] The following domain diameters will achieve such 
a limit for the associated LIq peak locations in the total 
structure factor (arranged by increasing domain diame- 
ter): 
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d = 27.5 A 
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d = 37.8 A 
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d = 53.4 A 


Qs + AQ+ 


d = 76.6 A 



(3) 

To illustrate the accuracy of such a simple analysis, 
consider the example of a NiMn crystal with a normal- 
ized distribution of spherical particle diameters d' (How- 
ell et aL 30 consider an entire family of distributions, but 
the characteristic peak broadening given by Eq. ^is in- 
dependent of the actual distribution), 



p ^-U4 



(4) 



In the limit of very large average domain diameters, 
the Bragg peaks from both structures are easily resolved 



with half the domains comprising an fee structure and 
half an Llo structure. The average diameter of the do- 
mains is d and the width of the distribution is a = d/2. 
In this approximation, the regions between the spherical 
domains is a highly strained structure with no long-range 
order, and the domains are entirely uncorrelated with 
each other. The structure has a constant local atomic 
density, approximately, within all the domains and the 
regions in between. 

It can be shown that the PDF of such a structure is 
the average of the PDF of the two uncorrelated domains, 
G(r) = l/2[G/cc(r) + G Llo (r)]f e (r,d), where the enve- 
lope function for the distribution of particle sizes given 
by Eq. His f e (r, d) = (1 + 2r/d) exp(-4r/d). Note that 
an ideal PDF constructed in this way still exhibits per- 
fectly resolved pair distances for < r < d, so that a 
local probe such as XAFS should readily be able to de- 
termine the presence of the two structures, and their rel- 
ative amounts. 

The structure factor is obtained from the PDF by a 
Sine transform, and is thus the average of the structure 
factors of the two uncorrelated sets of domains. The lat- 
ter can be obtained by convoluting the Cosine transform 



of the PDF envelope function with the structure factors 
of the associated infinite crystals. The expression for the 



total structure factor is then 
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where 
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is the convolution function for the distribution given by 
Eq. 21 It has a Lorentzian-like shape that serves to 
broaden the Bragg peaks in the structure factors of the 
ideal structures. 
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FIG. 4: The structure factor S(Q) for a NiMn alloy with equal 
amounts of Llo and fee structures within embedded domains 
with the average diameters shown. The vertical solid line is 
the first ideal Bragg peak position shared by the two coex- 
isting structures. The two vertical dashed lines are the addi- 
tional ideal peak positions for the fee structure, and the four 
vertical dotted lines are the additional ideal peak positions 
for the Llo structure. The splitting widths, given by Eq. [5] 
are labeled at the top. 

Figure 0] shows the structure factor of the NiMn alloy 
for each of the average diameters given in Eq. obtained 
by convoluting the ideal structure factor (seven delta 
functions representing the ideal Bragg peaks of the fee 
and LIq phases) with the function given by Eq. [fJl Also 
shown is the case when d — 100 A, where all peaks are 



easily resolved. The vertical lines are the locations of the 
ideal Bragg peak positions for the Llo and fee structures, 
and serve as guides to determine when individual peaks 
are resolvable. Indeed, for d — 76.6 A, 53.4 A, 37.8 A, 
and 27.5 A, the seventh, second, fifth, and fourth peaks 
become entirely unresolved, respectively. For all average 
domain sizes d < 30.0 A, the total structure would be 
characterized as entirely fee if one considered only the 
positions of peak maxima and shoulders in S(Q), while 
a local probe would readily resolve the two structures. 

The diffraction limit for uncorrelated domains can be 
determined with little analytical effort, and the domain 
size by which this limit is attained is quite small. When 
domains coexist but remain correlated, such that collec- 
tions of atoms that span across domain interfaces main- 
tain the coherence of the crystal, then the stresses intro- 
duced by differing orientations of chemical ordering can 
introduce significant strain throughout the domains, and 
such analysis is no longer feasible. For this scenario, we 
consider a numerical treatment to investigate the rela- 
tionship between domain size and diffraction limit. 



IV. THE STRUCTURE FACTOR AND PAIR 

DISTRIBUTION FUNCTION OF TWO 

DIMENSIONAL CORRELATED FCC AND Ll 

RECTANGULAR DOMAINS 



The two dimensional lattices represent [111] planes of 
the NiMn alloy films discussed, within which the Llo 
structure is comprised of alternating rows of the two el- 
ements. In one class of simulations, rectangular domains 
of various sizes of the LIq structure are embedded in a 
chemically disordered triangular (fee [111]) matrix, with 
each structure present in equal amounts. Rectangles are 
constructed with N rows and N columns of atoms, giv- 
ing an average side of length d = N(2.63 A+\/3/2 • 2.56 
A)/2. Another class of simulations considers hexagonally 
arranged domains of Llo chemical ordering that populate 
the entire crystal (no disordered fee component), with 
a diameter d that measures the span between opposing 
hexagonal vertices. 

Atomic interactions described with a Rose potential^ 
determine the precise atom positions required to mini- 
mize the total elastic energy at zero pressure. The equi- 
librium bond lengths are inferred from the observed three 
dimensional fee and tetragonal lattice parameters de- 



tailed in the last section: a = Ni-Ni = Mn-Mn = 2.63 A 
and b = Ni-Mn = 2.56 A, and a cutoff distance of 3.60 A 
limits the interactions to nearest neighbors. The chemi- 
cally disordered phase averages to an fee [111] structure 
with average bond length of a — 2.594 A. The bond ener- 
gies and moduli are equal for all three interactions, thus 
simplifying the model significantly while still capturing 
the relevant physics. The lattices are constructed with 
approximately one million atoms, and periodic bound- 
ary conditions are used in all calculations. Large lattice 
sizes are required to obtain the desired resolution in the 
diffraction data (Sq = 0.002 A -1 ), and as such three- 
dimensional simulations have exceeded our current com- 
putational limits (although preliminary results indicate 
that our conclusions will remain qualitatively the same) . 
Finally, our results given here are averages obtained from 
ten lattices, each differing in their random arrangement 
of domains, but statistical differences in their structure 
factors are small. The ability to obtain a minimum elastic 
energy within such large lattices allows us to demonstrate 
the importance of the elastic response of the lattice. 



The lattices are distinguished by the domain sizes of 
the Llo structure. Regions of typical lattices are depicted 
in Fig. |21 For the rectangular domains, portions of three 
domains show the possible orientations of the LIq [111] 
chemical ordering as well as the chemically disordered 
fee structure that separates them. We compare the cases 
in which no dislocations can form at the interfaces be- 
tween these phases, by maintaining the interactions of 
each pair of atoms that is determined by the perfect fee 
arrangement of the atoms prior to the energy minimiza- 
tion routine, to those in which dislocations are allowed by 
continually allowing for changes to the pairing of atomic 
interactions based on new nearest neighbor arrangements 
that occur during the crystal relaxations. The hexagonal 
domains represent the extreme case in which the crystal 
is entirely chemically ordered, and no dislocations can 
occur at any interfaces, such that the preferred tetrag- 
onal configurations within each domain present signifi- 
cant stresses at these interfaces. One therefore expects 
this case to represent the limit to which a large degree 
of related strain can occur throughout the crystal. In 
all cases, because the different domain orientations are 
present in equal amounts, the material does not exhibit 
any gross anisotropic distortion, and thus any changes 
in the average crystal length along both directions are 
equal. 



In the interval < Q < 5 A 1 , the ideal triangular 
structure in the [111] plane of the NiMn alloy has three 
Bragg peak positions at Qi = 47r/v3a, Q 2 — 47r/a, and 
Q3 = 8-K/VSa. The ideal Ll structure has seven Bragg 
peaks at Q 1 , Q 2 - AQ 2 , Q 2 , Q 2 + AQ 2 , Q 3 - AQ^ , 
Q3, and Q3 + AQ3 , where the splitting widths from the 



second and third peaks are given by 



AQ 2 " = 



AQ2 



47T 

7s ! 


(S- 


V a 2 ab b 2 


4-7T 


m- 


ab b 2 a 



0.0661 A 1 



0.0646 A l 



.0L(l-l\. „,„ 770 a- 1 



"* - 3s (s-i) -■""»*"'■ (7) 

Assuming that Eq. ^ might apply in two dimensions 
as well, we can at once conclude that if the two co- 
existing structures were entirely uncorrelated with each 
other, then when the average domain side lengths are 
larger than d = 3.48/AQj = 53.9 A, all of the struc- 
ture factor peaks would be easily resolved. When they 
become smaller than d = 3A8/AQ^ — 45.2 A, all peaks 
from the LIq [111] structure would be unresolved and the 
structure factor would indicate that only an fee [111] was 
present. Note the exceedingly small domain size interval 
between these two limits, which is only a consequence of 
considering two dimensional structures. 

Correlated structures preclude the use of such simple 
analysis, however, as the crystal structure within the Llo 
chemically ordered domains do not exhibit the tetragonal 
distortions seen in their uncorrelated counterparts. Re- 
quiring that every row in the two dimensional structure 
be coherent across the entire lattice introduces a signifi- 
cant energy constraint into the system. We shall see that 
the behavior of the resulting strain is somewhat nonintu- 
itive, by consistently changing the ordered Ll structure 
into something very different than that suggested by the 
pair interaction potentials alone. Indeed, considering the 
two cases (50% and 100% Llo ordering) offers two ex- 
treme consequences of such strain. 

Figure[S]shows the structure factors of the NiMn alloys 
for the correlated rectangular domains (bottom) with 
sides of length d = 560 A, 280 A, and 140 A, includ- 
ing the cases in which dislocations are allowed. In all 
correlated cases (no dislocations), the second, fourth, 
and fifth Bragg peaks for the tetragonal structure are 
never present (the vertical lines in the figure represent 
the ideal Bragg peak positions of both structures). Con- 
sidering the analysis of the three dimensional structures 
in the previous section, it is clear that their absence is 
not strictly a matter of the structure having attained 
the diffraction limit. Instead, the tetragonal structure 
is not able to form entirely within the domains due to 
the coherence condition at their interfaces, and thus the 
structure factor is revealing an average structure of the 
fee and LIq structures. When dislocations are present, 
the tetragonal structure is easily resolved at these do- 
main sizes, and the peaks show a trend consistent with 
their three dimensional counterparts, in which decreasing 
domain sizes diminishes the peak heights at the expense 
of increased peak width. One cannot conclude, however, 
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FIG. 5: The structure factor S(Q), in arbitrary units and off- 
set for comparison, for two dimensional [111] planes of NiMn 
alloy. Bottom: the structure factor from lattices with equal 
amounts of fee and Llo structures, the latter in embedded 
rectangular domains with the average side lengths shown. 
Top: the structure factor from lattices entirely populated with 
hexagonal Llo ordered domains, scaled by a factor of seven 
for comparison. The three vertical dashed lines are the ideal 
Bragg peak positions shared by fee and Llo structures. The 
four vertical dotted lines are the additional ideal peak posi- 
tions for the Llo tetragonal structure. The splitting widths 
are labeled at the top. 



if any fee structure is present in such systems, as the fee 
Bragg peaks coincide with those of the tetragonal struc- 
ture (Qi, Qi, and Qz) by construction. A real-space 
treatment, utilizing the pair distribution function G(r) 
will address this issue. 

The upper portion of Fig. [3] shows S(Q) for the com- 
plete chemically ordered crystals with correlated hexago- 
nal domains. The peak structure is consistent with that 
of a single mesoscale glassy phase, with extreme peak 
broadening due to short length-scale structural ordering. 
In addition, significant shifts of the peak positions from 
the tetragonal phase, compared to structures with rect- 
angular domains, indicate the presence of a large tetrag- 
onal distortion. Indeed, the hexagons in Fig. [21 are ori- 
ented in a consistent manner such that no two neighbor- 
ing domains share the same orientation, which promotes 



tetragonal distortions far larger than those dictated by 
the interaction potentials. This is particularly noticeable 
from the shifts in the peak positions about Qi and Q$. 
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FIG. 6: The pair distribution function G(r), in arbitrary units 
and offset for comparison, for the same lattices represented in 
Fig. The three vertical dashed lines denote are the ideal 
pair distances for the fee structure, while the six vertical dot- 
ted lines are for pair distances of the tetragonal structure. 
G(r) from lattices entirely populated with hexagonal Llo or- 
dered domains are scaled by a factor of two for comparison. 

To better understand the structures, we show in Fig. 
the pair distribution function G(r) for the same relaxed 
lattices, with the dashed vertical lines denoting the ideal 
pair distances for the fee structure and the vertical dotted 
lines for those of the tetragonal structure. Interestingly, 
when dislocations are allowed at the rectangular inter- 
faces (bottom), the entire lattice seems to relax to the 
tetragonal structure, even though the Llo chemical or- 
dering exists in only half the crystal. Only a very low 
amplitude peak is seen at the third fee pair distance in 
these cases. This may be a result of collective behavior 
among the atoms within the domains, in which the order- 
ing allows for a common, global interaction-like scenario 
that easily influences the random interactions outside the 
domains. When coherence across the interfaces is en- 
forced, however, the distribution of lengths from the first 
through third nearest neighbors generally indicates the 
presence of only one average fee structure, with tetrag- 
onal peak splittings being revealed only in the smallest 
domain size considered. 



The lattices with entire hexagonal domains exhibit ex- 
tremely broad and shifted peaks in the pair distribution 
function, a result of significant and continuous tetrago- 
nal distortions. The peak structure present in G(r) de- 
cays away to a uniform atomic density at approximately 
r = 25 A (not shown in Fig. EJ), which suggests a descrip- 
tion of the structure as a mesoscale glass with medium- 
range order. Additionally, the peak positions, when re- 
solvable, are shifted from the LIq tetragonal structure so 
much that one can only conclude that some collective be- 
havior among the domains is present, which we alluded 
to earlier with respect to their structure factors. 

When there is only one average structure present 
throughout the crystal, such as the case with correlated 
rectangular domains of LIq chemical ordering compris- 
ing half the structure, then local and extended probes 
reveal the same structure. Then, in the context of the 
NiMn alloy considered here, magnetic order is probably 
more correlated with the LT chemical ordering than the 
tetragonal distortions. 



V. CONCLUSION 

We have provided conclusive evidence, both analyti- 
cally and numerically, that it is indeed possible to have 
well-defined tetragonal distortions within localized do- 
mains that fall below the diffraction limit. The structure 
is therefore readily attainable via local probes such as 
XAFS. An analytical treatment considered the phase co- 
existence of such structures with a chemically disordered 
fee structure, where the two structures are entirely un- 
corrected. Basic diffraction analysis predicts with great 



accuracy the domain size below which the tetragonal sig- 
natures are no longer present in unique Bragg peaks, but 
instead become part of the diffuse scattering of the aver- 
age fee structure. 

In the case in which atomic coherence is maintained 
across the interfaces of domains, numerical treatments 
incorporating elastic energy considerations indicate that 
the tetragonal structure may only manifest itself far away 
from the interfaces, but the domain sizes considered in 
this treatment never allowed for true fee and tetragonal 
phase separation. Instead, the strain incurred by the 
competing domains is of such great extent as to produce 
even larger tetragonal distortions than the pair interac- 
tion potentials dictate, as well as significant disorder in- 
dicative of mesoscale glassy materials (seen in the 100% 
LIq ordered hexagonal domains). Alternatively, it has a 
tendency to average the structure to an fee phase (as in 
the case of 50% LTo ordered rectangular domains). 

If the latter case is indeed the correct description of the 
nanoscale ordering of NiMn alloys, then one must con- 
clude that the presence of antiferromagnetism is a direct 
result of the chemical ordering present in the LIq struc- 
ture, and not the actual tetragonal distortion since this is 
suppressed. This type of smooth transition in the bond 
lengths is however what was observed experimentally in 
the NiMn system^ 
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